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ABSTRACT 

Results  are  presented  from  the  study  of  the  off-  axis  radiation  field 
caused  by  a  relativistic,  mono-directional  electron  beam  passing  through 
water.  Using  the  Naval  Postgraduate  School  100  MeV  linear  accelerator, 
off-axis  radiation  dose  was  measured  with  calcium  fluoride  thermo- 
luminescent dosimeters  placed  at  various  path  lengths  out  to  two  radiation 
lengths.  The  off-axis  dose  was  calculated  using  the  electron  transport  code 
CYLTRAN  of  the  integrated  TIGER  series  of  coupled  electron/photon  Monte 
Carlo  transport  codes.  Calculations  were  performed  at  Los  Alamos  National 
Laboratory.  Comparison  with  the  results  is  made  and  CYLTRAN  is  found  to  be 
in  agreement  with  experimentally  measured  values.  The  extension  of 
experimental  results  in  one  medium  (water)  to  another  (air)  appears  to  be 

valid. 
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I.      INTRODUCTION 

Recent  development  of  new  high  energy,  high  current  accelerators  has 
renewed  interest  in  the  ability  to  accurately  predict  radiation  exposure  in  the 
vicinity  of  a  monodirectional  electron  beam.  In  addition  to  personnel  safety 
requirements,  the  areas  of  vulnerability  and  lethality  as  they  apply  to 
charged  particle  beams  are  becoming  increasingly  important.  The  majority 
of  the  literature  on  radiation  dosimetry  applies  to  low  or  moderate  energy 
sources.  The  experimental  verification  of  the  calculations  as  they  apply  to 
high  energy  electron  beams  has  been  sparse  and  comparison  of  several 
computational  methods  has  been  inconclusive. 

Recent  experiments  at  the  Naval  Postgraduate  School  by  P.F.  Cromar 
attempted  to  address  some  of  these  questions  by  studying  the  off-axis 
radiation  field  in  liquid  nitrogen  as  a  function  of  radiation  length  [Ref.  I:  p. 4]. 
Focused  electrons  of  100  MeV  incident  energy  were  delivered  to  insulated 
containers  filled  with  liquid  nitrogen.  The  measured  dose  was  compared  to 
calculations  using  the  electron  transport  code  CYLTRAN  .  Cromar's  results 
compared  favorably  (within  a  factor  of  2)  with  the  results  predicted  by  the 
calculations.  However,  where  comparisons  could  be  made,  Cromar's 
results  appeared  to  differ  from  an  earlier  calculation  by  R.A.  Lindgren 
using  the  code  ETRAN-16, which  referred  to  50  Mev  electrons  in  air 
[Ref.  I:  p. 39]  .  It  is  the  objective  of  this  thesis  to  address,  through  a  series 
of  experiments,  several  questions  concerning  the  validity  of  calculations  and 


experimental  methods  as  they  apply  to  high  energy  electron  beams.    Passive 
dosimetry  will  be  used  to  determine: 


1)  Can  the  computer  computations  accurately  predict  the  experimental 
results? 

2)  Can  the  results  in  one  medium  be  used  to  predict  the  response  in  a 
second  medium? 

3)  Is  there  an  energy  dependence  to  the  response? 

4)  If  the  responses  are  different,  what  is  the  significance  as  it  applies  to 
the  energy  spectrum? 


II.      BACKGROUND 

A.    INTERACTION  OF  RADIATION  WITH  MATTER 

When  a  primary  electron  passes  through  an  absorber  it  will  only  have  a 
well-defined  range  for  energies  of  less  than  a  few  MeV.  At  these  low 
energies,  energy  loss  is  by  ionization  and  excitation  of  the  atomic  electrons 
in  the  atoms  of  the  absorber[Ref.2:  p. 508].  The  energy  loss  of  electrons  due 
to  ionization  is  given  by  [Ref.3:  p. 38]: 


'dp,  =    4-rte4  n 
•ljx'ion    mc2 


?mr2      3  1  %      1  1 
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where     e  =  charge  on  the  electron    (cgs) 
n  =  electron  density 
m  =  mass  of  electron 
c  =  speed  of  light 
I  =  ionization  potential 
6  =  ratio  of  electron  velocity  to  c 

As  the  electron  energy  increases,   radiative  energy  loss  becomes  more 
important  than  ionization.   This  energy  loss  is  given  by  [Ref.  3:  p. 68]: 
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where     r0  =  classical  radius  of  the  electron 

h  =  Planck's  constant 
Z  =  atomic  number 

u   =    frequency   of   the   highest   energy 
photon  produced 

The  energy  loss  by  radiation,  called  bremsstrahlung,  is  proportional  to 

Z2  of  the  material  and  increases  linearly  with  the  electrons  energy,  whereas 

losses  due  to  ionization  and  excitation  are  only  proportional  to  Z.     Thus  the 

radiation  loss  predominates  at  the  higher  energies.      The  approximate  ratio 

of  the  two  losses  is  [Ref.3:  p. 64] 
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There  is  a  critical  energy,  Ec,  for  which  the  two  energy  losses  are  equal: 


800 
E,-.  (MeV)  =  —t- 


If,  in  a  given  material,  the  incident  energy  is  much  larger  than  the 
critical  energy,  then  radiation  loss  predominates  over  ionization  loss  and  can 
be  considered  to  be  the  only  significant  loss  mechanism  [Ref.3:  p. 72]. 


When  electrons  of  high  initial  energy  pass  through  matter,  they  lose 
their  energy  mainly  by  emitting  bremsstrahlung  photons.  These  high-energy 
photons,  or  gamma  rays,  in  turn  lose  their  energy  mainly  by  the  electron 
pair-production  process.  By  a  series  of  similar  interactions  and  a 
succession  of  these  simple  elementary  processes,  a  multiplication  process 
comes  about.  An  avalanche  of  electrons  and  photons  develops  and  is  referred 
to  as  cascade  shower  or  electron-photon  cascade. [Ref. 4:  p. 91  ][Ref . 5:  p. 221] 

The  energy  of  the  primary  electron  is  divided  among  the  cascade 
particles,  their  sum  equaling  the  primary  energy  except  for  the  amount  lost 
due  to  ionization.  A  rapid  multiplication  of  shower  particles  results  with  the 
number  of  particles  growing  exponetially  as  long  as  their  energies  exceed  the 

critical  energy,  Ec.      As  the  number  of  particles  reaches  a  maximum  value 

and  the  energy  per  particle  is  smaller,  the  ionization  losses  become  more 
important    than    the    bremsstrahlung    losses.  Multiplication    of    shower 

particles  ceases  and  as  energy  is  spent  in  ionization  and  excitation  processes 
the  number  of  particles  decreases.  [Ref. 4:pp. 92-93] 

The  maximum  number  of  shower  particles  amounts  to  : 


■''max     Ec 


This  maximum  is  observed  at  an  absorber  depth  of: 
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t  ln(Eo/Ec) 

'■max  =  

In  2 

with  tmax  being  measured  in  radiation  lengths. [Ref. 4:  p. 93] 

Photon  production  by  bremsstrahlung  and  their  angular  distribution  is 
largely  a  statistical  process.  The  spectrum  of  energies  produced  by  an 
electron  beam  will  be  broad  and  the  percentage  of  energy  transferred  to  any 
given  photon  following  an  interaction  will  follow  a  probability  distribution 
function  determined  by  the  energy  of  the  incident  electron  and  the  atomic 
number  of  the  material.  The  result  is  that  many  secondary  electrons  which 
are  emitted  cause  further  ionization  and  excitation,  and  may  penetrate  far 
beyond  the  range  of  the  primary  electrons  [Ref. 2:  p. 509].  The  precise 
mathematical  development  of  showers  is  complicated  and  can  contain  large 
fluctuations  from  the  average  behavior.  In  practice,  it  becomes  necessary 
to  use  Nonte  Carlo  techniques  to  calculate  the  shower  induced  dose. 

The  difficulty  in  separating  the  effects  of  the  primary  and  secondary 
electrons   makes   range     an  inappropriate  measure  of  electron  absorption 

above  a  few  MeV  energy  [Ref. 2:  p. 509.].  The  concept  of  radiation  length,  XQ, 

is  useful  in  the  interpretation  of  this  energy-loss  phenomena  and  is  defined  by 
Fenyves  and  Haiman  [Ref  .4:  p. 25]  as*. 


=   4c<|  Z(Z+1  )rn2[ln(183  Z      ")] 


1+n. 12    — 


Radiation  length  can  be  visualized  as  the  distance  in  which  the  electron 
energy    is  reduced  to  1/e  of  its  initial  value  and    is  independent  of  electron 
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energy.  In  practice  it  is  more  convenient  to  express  the  quantity  in  units  of 
gm/cm^  (Table  1)  in  which  the  calculated  values  form  a  smooth  curve  as  a 
function  of  Z.  [Ref.2:  pp. 510-512] 

TABLE  1 
VALUES  OF  RADIATION  LENGTH  FOR  VARIOUS  SUBSTANCES 


Substance 

Z 

A 

Radiation  lengths 
(gm/cm2)     (cm) 

Crit 

cal  energy 
(Mev) 

Carbon 

6 

12 

44.6 

30.0 

102 

Nitrogen 

7 

14 

39.4 

88.7 

Air 

7.37 

14. 

78 

37.7 

31. Ox 

103 

84.2 

Water 

7.23 

14. 

3 

37.1 

37.1 

83.8 

Oxygen 

8 

16 

35.3 

77.7 

12 


III.     ISSUES 

Previous  experiments  dealing  with  the  experimental  verification  of 
calculations  as  they  apply  to  high  energy  particle  beams  raised  several 
issues.  Direct  comparisons  between  experiments  and  calcualtions  have  been 
difficult  due  to  differing  experimental  geometries,  varying  incident  energies, 
and  multiple  computational  computer  codes.  This  experiment  will  address 
some  of  these  issues. 

A.    EXPERIMENTAL  GEOMETRY 

The  recent  experiments  by  P.F.  Cromar  [Ref.  I]  were  conducted  in  an 
insulated  container  filled  with  liquid  nitrogen.  The  use  of  liquid  nitrogen  as 
the  test  medium  precluded  the  immersion  of  the  dosimeters  in  the 
environmental  matrix.  The  dosimeters  were  placed  in  air  on  the  exterior  of  a 
one  inch  thick  hard  foam  container.  Ideally,  direct  comparison  with 
calculation  should  be  done  with  the  dosimeters  placed  inside  the  test  medium 
so  that  the  side  and  back  scatter  contributions  are  the  same  for  experiment 
and  calculations. 

Separate  tanks  were  required  for  measurements  at  different  distances 
from  the  beam  entrance.  As  a  result,  several  independent  runs  were 
necessary  to  compile  the  data  at  all  distances.  If  the  dosimeters  could  be 
placed    directly    in    the    test    medium    all    measurements    could    be    done 
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simultaneously  to  minimize  experimental  beam  condition  variability.  For  this 
experiment,  thermoluminescent  dosimeters  were  sealed  in  plastic  and  placed 
directly  in  water. 

B.  ENERGY  DEPENDENCE 

The  experiment  by  Cromar  was  conducted  at  a  single  energy  of  100  MeV. 
This  is  roughly  the  critical  energy  of  liquid  nitrogen  where  there  is  equal 
contribution  from  bremsstrauhlung  and  ionization  to  the  energy  loss 
mechanism.  The  contributions  to  energy  deposition  are  different  from  that 
for  energy  loss.  It  is  possible  that  the  continuous  slowing  down 
approximation  in  the  computational  code  may  not  be  adequate  at  all  energies. 
Experiments  were  conducted  at  100  MeV  and  20  MeV.  The  experimental 
results  from  the  20  MeV  run  were  not  returned  for  inclusion  is  this  paper,  but 
results  of  other  published  results  at  varying  energies  are  discussed. 

C.  COMPUTATION 

It  is  not  clear  from  the  published  data  and  the  calculations  of  Cromar  that 
all  versions  of  the  electron  transport  codes  produce  the  same  computational 
results.  The  apparent  disagreement  between  ETRAN-16  and  CYLTRAN  needs  to 
be  investigated.  In  translating  the  calculated  dose  in  liquid  nitrogen  to  that 
which  might  be  observed  in  air,  a  single  geometrical  extrapolation  was  used. 
This  may  not  be  adequate  as  multip^  scattering  angles  and  radiation  lengths 
in  the  beam  axis  and  the  perpendicular  direction  interplay. 
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In  this  experiment,  the  results  of  the  two-dimensional  CYLTRAN  code  for 
the  energy  deposition  in  water  have  been  converted  to  a  one-dimensional 
geometry  and  compared  with  the  published  results  for  ETRAN.  In  addition, 
CYLTRAN  results  for  water  and  liquid  nitrogen  are  scaled  and  compared  to 
the  CYLTRAN  predictions  for  air. 

D.     DOSIMETRY 

CaF2  dosimeters  have  a  non-linear  response  as  a  function  of  photon 
energy  below  200  keV  as  shown  in  Figure  3.  1 . 
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Figure  3.  1    Calculated  Response  of  CaF->  and  LiF 


Should  the  radiation  field  in  the  vicinity  of  the  electron  beam  result  in  a 
large  fraction  of  the  deposited  energy  being  due  to  lower  energy  photons, 
then  the  dose  measured  could  be  in  serious  error.    The  dosimeters  have  been 


calibrated  with  respect  to  a  ^°Co  source.  The  spectrum  due  to  a  line  electron 
source  is  not  identical  to  that  from  6oCo. 

Dosimeters  of  CaF2  and  LiF  were  placed  in  symmetrical  positons  for  a 

direct  comparison  of  the  measured  dose.  This  portion  of  the  experiment 
remains  incomplete  as  the  results  for  the  LiF  dosimeters  were  not  returned 
due  to  a  faulty  set  of  calibration  curves  at  the  reading  facility. 
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IV.     CALCULATIONS 

A.     DESCRIPTION 

The  calculations  of  electron/photon  showers  were  performed  using  the 
computer  program  CYLTRAN  of  the  Integrated  TIGER  Series  of  Transport 
Codes[Ref.  6].  The  purpose  of  the  calculation  was  to  determine  the  off-axis 
dose  resulting  from  discharging  a  monodirectional  beam  of  electrons  into 
several  different  media. 

The  cylindrical  geometry  of  the  computer  code  CYLTRAN  is  well  suited  for 
the  desired  output.  A  0.8  cm  diameter  electron  beam  of  20,000  particles 
was  incident  normal  to  targets  of  water,  air,  and  liquid  nitrogen.  Each 
target,  as  indicated  in  Figure  4.1,  was  divided  into  zones  with  energy 
deposition  being  recorded  at  distances  that  corresponded  to  0.0,  0.25,  0.50, 
1.0,  1.5,  and  2.0  radiation  lengths  in  each  medium.  At  each  distance,  ten 
concentric  sub-zones,  Figure  4.2,  were  created  to  measure  the  off-axis 
energy  deposition.  Zone  thickness  was  kept  thin  to  provide  semi-discrete 
data  points  at  specific  distances  of  penetration  into  the  medium.  However, 
sufficient  thickness  was  necessary  to  allow  for  the  finite  track  segments  used 
in  the  computer  calculations  and  to  provide  for  statistical  validity  of  the 
results. 

The  output  of  the  code  included  the  target  material,  zone  or  sub-zone 
location  and  mass,  energy  deposited  by  primary  and  knock-on  electrons, 
energy  deposited  by  photons,    and  the  total    energy  deposited  within   each 
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Figure  4.1     CYLTRAN  Geometry 
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Figure  4.2  CYLTRAN  Sub-Zones  for  Water 


sub-zone.  The  simulation  was  run  in  ten  batches  of  2000  particles  each  with 
statistical  values  being  computed  after  each  batch.  Each  value  of  energy 
deposition  is  followed  by  a  two  digit  number  that  represents  the  one-sigma 
uncertainty  that  results  from  the  Monte-Carlo  routine,  with  low  numbers 
indicating  narrow  bounds  of  uncertainty. 

At  distances  well  within  the  targets  the  statistical  uncertainties  grew  due 
to  the  fewer  number  of  particles  with  sufficient  energy  to  penetrate  to  these 
zones.  Decreasing  the  statistical  uncertainties  can  be  accomplished  by 
increasing  the  number  of  incident  particles,  however  the  relationship 
between    numbers    of    particles    and    uncertainties    is    not    linear.  Each 

simulation  run  of  20,000  particles  used  about  twenty  minutes  of  Cental 
Processor  Unit  time  on  the  Cray  computers  at  Los  Alamos.  Increasing  the 
number  of  particles  to  45,000  would  have  reduced  the  uncertainty  by  only 
about  ten  percent  but  would  have  increased  the. CPU  time  by  a  factor  of  two. 

Since  the  primary  purpose  of  the  computer  simulations  was  to  provide 
immediate  comparisons  with  the  experimental  results,  no  attempt  was  made 
to  fully  optimize  the  numerical  calculations.  When  reasonable,  all  default 
values  within  the  codes  were  used.  This  led  to  track  segments  that  were 
larger  than  ideal  at  the  higher  energies  near  the  entrance,  while  the  track 
segments  were  too  small  for  the  less  energetic  particles  well  within  the 
target.  It  is  believed  that  for  the  comparisons  desired,  the  integration  of 
the  large  number  of  incident  particles  is  adequate  to  provide  meaningful 
results  inspite  of  the  less  than  optimum  analysis. 
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Appendix  A  contains  a  description  of  CYLTRAN  and  sample  input  data 
streams  for  the  generation  of  the  cross-sections  (Table  A1)  and  the 
execution  of  the  program  (Table  A2).  The  CYLTRAN  output  of  the  energy 
deposition  normalized  to  one  incident  particle  for  100  MeV  electrons  in  water, 
air,  liquid  nitrogen,  and  20  MeV  electrons  in  water  is  contained  in  Tables 
A3-A25. 

B.     RESULTS 

The  total  energy  deposited  in  each  sub-zone  was  converted  to 
Rads/Coulomb  for  comparison  with  the  experimental  results.  The  normalized 
dose  due  to  100  MeV  electrons  for  each  sub-zone  in  water  is  contained  in 
Table  B1  of  Appendix  B.  The  resultant  graphs  of  the  off-axis  dose  for  each 
penetration  depth  are  given  in  Figures  Bl  through  B7.  Appendix's  C,  D,  and  E 
contain  similar  results  for  the  100  MeV  electrons  in  air  and  liquid  nitrogen, 
and  the  20  MeV  electrons  in  water. 

The  curves  drawn  on  the  graphs  are  not  intended  to  necessarily  be  the 
mathematically  best  fit  approximation  for  the  data  although  in  most  cases  the 
difference  is  minimal.  Instead,  the  curves  were  generated  using  the 
exponential  and  polynomial  equations  of  Appendix  F,  in  order  to  construct  the 
three  dimensional  half-axis  surface  plots  of  Figures  4.3  through  4.6  and  the 
contour,  or  isodose,  curves  of  Figures  4.7  through  4.10.  The  vertical 
error  bars  result  from  the  indicated  one-sigma  errors  generated  by 
CYLTRAN.       The  plots  have  been  smoothed  only  in  the  off-axis  direction. 
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Smoothing  parallel  to  the  beam  axis  would  have  involved  a  good  deal  of 
numerical  experimentation  and  iteration  and  would  not  have  added  any 
meaingful  interpretation  to  the  results.  In  spite  of  the  modest  smoothing, 
the  relatively  small  number  of  data  points  resulted  in  the  rough  surface  plots 
of  Figures  4.7  through  4. 10  and  it  is  felt  that  no  additional  smoothing  could  be 
justified  based  on  these  few  points. 
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V.     EXPERIMENTS 

A.     DISCUSSION 

All  experimental  measurements  were  conducted  at  the  Naval 
Postgraduate  School  Accelerator  Lab  using  the  100  MeV  linear  electron 
accelerator.  A  rectangular  container  (100  x  46  x  38  cm3  )  of  4  mm 
polyethylene  plastic  was  used  to  contain  the  experimental  matrix.  This 
allowed  for  a  usable  test  area  of  10  cm  on  either  side  of  the  central  axis  with 
a  minimum  of  9  cm  of  water  beyond  this  to  provide  a  uniform  scattering 
medium  without  side  effects.  The  length  of  the  tank  allowed  measurements 
out  to  two  radiation  lengths  in  water. 

Calcium  fluoride  thermoluminescent  dosimeters  (TLDs)  were  provided  by 
Naval  Surface  Weapons  Center.-  The  TLDs  were  wrapped  in  a  single  layer  of  I 
mil  aluminum  and  were  enclosed  in  a  thin  plastic  film  to  prevent  moisture 
from  contaminating  the  crystals.  The  TLDs  were  mounted  to  soft  wood 
stretchers  at  intervals  indicated  in  Figure  5. 1  and  Figure  5.2.  Wood  contains 
hydrogen  and  carbon  and  is  similar  to  water  in  atomic  properties.  Pairs  of 
TLDs  (one  CaF2  and  one  LiF )  were  mounted  in  symmetric  positions.  The  TLDs 
were  generally  located  at  one  centimeter  spacing,  but  alternating  between 
sides  of  the  central  axis,  i.e.  1,3,5  cm  right  and  2,4,6  cm  left  .  This 
served  to  reduce  interference  between  adjacent  dosimeters  and  was  intended 
to  provide  a  means  of  interpreting  the  results  in  case  of  a  slightly  off-axis 
beam  or  drifting  of  the  stretcher  from  centerline.     TLD  positions  were  known 
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Figure  5.  1     TLD  Positioning  Within  Test  Tank  for  100  Mev  Run 
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to  within  0.5  cm  on  the  stretcher  while  the  stretcher  position  was  known  to 
within  one  centimeter. 

Pulsed  beams  of  100  MeV  and  20  MeV  electrons  were  normally  incident  on 
the  target  at  a  repetition  rate  of  60  Hz.  The  beam  passed  through  a  vacuum 
test  chamber  which  was  sealed  with  a  five  mil  aluminum  window  and  then 
through  25  cm  of  air  prior  to  entering  the  test  tank. 

The  total  charge  delivered  was  monitored  by  a  Secondary  Emission 
Monitor  (SEM).  A  vibrating  reed  electrometer  coupled  with  a  one 
micro-farad  capacitor  was  used  to  determine  the  charge  delivered  to  the 
SEM.  With  an  SEM  efficiency  of  .028,  the  total  charge  delivered  to  the  matrix 
could  be  determined.  Due  to  uncertainties  in  this  efficiency  factor,  it  is 
estimated  the  total  charge  delivered  can  only  be  determined  to  within  five 
percent.  Electron  energy  was  controlled  by  the  use  of  a  magnetic  deflection 
system  and  is  considered  to  be  sufficiently  stable  that  any  fluctuations  would 
not  be  noticable  within  the  limits  of  this  experiment. 

Phosphor  screens  were  placed  at  both  ends  of  the  tank  on  the  central  axis 
to  align  the  beam.  The  TLDs  were  removed  and  the  beam  was  passed  through 
the  empty  tank.  By  remotely  observing  the  phosphor  screens,  an  attempt 
was  made  to  make  the  beam  as  parallel  as  possible  while  within  the  tank. 
The  beam  diameter  entering  the  tank  averaged  0.5  cm  as  indicated  by  the 
phosphor  screens,  while  beam  divergence  caused  the  exiting  diameter  to 
approach  2.0  cm. 
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Once  alignment  and  focussing  were  completed  the  phosphor  screens  were 
removed,  TLDs  Installed,  and  the  tank  filled  with  water.  To  prevent 
saturation  of  the  TLDs  close  to  the  beam  axis  and  yet  still  provide  a 
meaningful  dose  on  the  distant  TLDs,  three  irradiation  segments  were 
conducted  for  each  run.  Based  on  estimated  dose  calculations,  selected 
TLDs  were  removed  after  each  segment  to  prevent  overexposure. 
Normalization  of  the  data  to  Rads  per  Coulomb  delivered  will  account  for  the 
different  exposure  times. 

B.    RESULTS 

The  TLDs  were  mailed  to  the  Naval  Surface  Weapons  Center,  White  Oak, 
Maryland  for  reading.  The  results  of  the  TLD  readings  were  reported  in  Rads 
and  using  the  recorded  charge  delivered  to  each  TLD,  the  readings  were 
normalized  to  Rads  per  Coulomb.  For  the  100  MeV  experiment  a  faulty  set  of 
calibration  curves  prevented  proper  reading  of  the  LiF  dosimeters  and  these 
results  were  not  returned. 

The  normalized  dose  from  the  CaF2  dosimeters  is  contained  in  Table  Gl  of 

Appendix  G.  Figures  5.3  through  5.8  are  plots  of  the  normalized  off-axis 
dose  at  each  length  within  the  water.  The  vertical  error  bars  are  drawn 
based  on  a  ten  percent  error  in  the  normalized  dose  .  The  horizontal  error 
bars  are  the  result  of  the  uncertainty  in  the  TLD  position  on  the  stretcher 
(±0.5  cm),  the  strecher  location  with  respect  to  the  central  axis  of  the  test 
tank(±  1 .0  cm),  and  the  skewness  of  the  beam  with  respect  the  center  line  of 
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the  tank  (±2.0  degrees).  The  curves  drawn  on  the  figures  are  the 
semi-smoothed  results  from  the  CYLTRAN  calculations  and  are  included  for 
comparison. 

The  objective  of  the  experiment  was  to  provide  quick  but  meaningful 
results  for  comparison  with  the  calculations.  As  a  result,  the  experimental 
apparatus  was  kept  simple  and  straight  forward.  As  a  consequence,  the 
cumulative  errors  at  the  extreme  distances  result  in  horizontal  error  bars 
that  encompass  a  large  portion  of  the  horizontal  axis.  However,  by  placing 
sequential  TLDs  on  alternate  sides  of  the  central  axis,  some  simple  yet 
justified  and  reasonable  corrections  can  be  applied  to  the  data  in  order  to 
present  the  true  picture. 

Figure  5.3  shows  how  the  true  off-axis  position  of  the  TLDs  varies  with 
the  degree  of  skewness  in  the  beam.  Since  alternate  alternate  TLDs  were 
placed  on  alternate  sides  of  the  matrix  x-axis,  this  skewness  results  in  a 
displacement  of  the  true  off-axis  distance  measured  from  the  centeriine  of 
the  beam. 

Table  G2  of  Appendix  G  contains  a  set  of  corrected  positions  for  the  TLDs 
based  on  varying  degrees  of  skewness  for  the  beam.  A  misalignment  of  only 
two  degrees  from  parallel  results  in  an  error  of  the  off-axis  position  of  2.58 
cm  at  the  furthest  distance  for  which  data  was  recorded.  Figures  5.9 
through  5.12  present  the  results  if  a  uniform  correction  of  1.7  degrees  is 
applied  to  all  the  data  points.  At  the  inner  distances  the  uncertainty  in  the 
position  of  the  stretcher  is  greater  than  the  correction  of  the  beam  skewness 
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Figure  5.3  :  Effect  of  Beam  Skewness  on  Off-axis  Zero  Point 
and  the  effect  of  applying  the  skew  correction  is  not  as  visable.    In  all 
cases,   the  data  generally  alternates   between  high   and  low  points   of  the 
calculated  curves,  in  relation  to  being  on  the  right  side  or  the  left  side  of  the 
central  axis. 
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VI.    CONCLUSIONS 

A.    COMPARISON  OF  COMPUTER  CALCULATIONS  AND  EXPERIMENTAL  RESULTS 

Within  the  limits  of  this  experiment,  there  is  good  agreement  between  the 
measured  and  calculated  normalized  dose  due  to  100  MeV  electrons  in  water. 
The  experimental  results  at  0.0  cm  and  9.2  cm  were  significantly  higher  than 
predicted.  The  TLDs  at  0.0  cm  were  intended  to  measure  the  dose  as  a  result 
of  backscattering.  However,  the  entrance  to  the  test  tank  was  not  shielded 
and  it  is  probable  that  these  front  TLDs  were  measuring  a  significant  amount 
of  radiation  scattered  from  the  aluminum  window  covering  the  exit  of  the 
vacuum  chamber. 

The  results  could  be  greatly  improved  with  a  few  simple  modifications  to 
the  TLD  matrix.  Alternating  consecutive  TLDs  from  one  side  of  the  central 
axis  to  the  other  showed  that  the  beam  was  not  directly  on  that  axis.  To  adjust 
for  this  discrepancy  it  was  necessary  to  try  several  iterative  corrections  of 
beam  angle  to  bring  the  data  back  in  line.  If  two  TLDs  had  been  placed  at 
symmetrical  positions  on  either  side  of  the  central  axis,  the  two  readings 
could  be  compared  ,  providing  a  means  of  calibrating  the  off-axis  positions  of 
the  TLDs  with  respect  to  the  true  beam  position.  In  addition,  use  of  a  rigid 
frame  to  support  the  TLDs  instead  of  the  floating  strechers  that  were  used, 
would  allow  for  much  more  accurate  positioning  of  the  entire  array. 
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B.     EXTENSIONS  OF  THE  PREDICTIONS  FROM  ONE  MEDIUM  TO  ANOTHER 

The  results  were  recorded  at  distances  which  correspond  to  fractions  of 
the  radiation  length  in  the  given  medium.  To  extend  the  results  from  one 
medium  to  another,  It  is  necessary  to  account  for  the  change  In  size  and  mass 
of  the  detector  as  the  radiation  length  varies  [Ref.1:  p. 34].  If  a  detector 
subtends  the  same  solid  angle  and  is  located  at  an  equivalent  distance  when 
measured  in  radiation  lengths,  the  same  energy  should  be  deposited  in  the 
detector. 
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Figure  6. 1    Ratio  of  the  areas  of  two  detectors 

As  shown  If  Figure  6. 1 ,  the  new  detector  has  an  area  and  hence,  a  mass, 
that  is  larger  by  a  factor  of  the  square  of  the  ratios  of  the  distance  to  the 
detector.     Since  dose  is  related  to  the  mass  of  the  detector, 

1  Rad  =  100  ergs  /  gm 
the  dose  at  the  new  detector  will  be  decreased  by  the  same  ratio.     Water  has 
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a  radiation  length  of  37.0  cm  and  air  has  a  radiation  length  of  307  meters,  so 
for  this  experiment  the  ratio  is: 


/R,  Y        /  37  \2  -6 

=  1 .45  x  1 0 


lR2,i    =    \ZQl) 


Table  6.1  contains  values  of  the  dose  in  air  predicted  from  the  dose  in 
water  compared  to  the  dose  in  air  predicted  by  CYLTRAN.  At  the  near 
distances,  the  results  vary  by  about  thirty  percent,  but  this  is  largely  due  to 
the  high  one-sigma  errors  of  the  data  at  these  distances.  At  one  radiation 
length  and  above,  the  statistical  errors  in  the  initial  data  are  much  smaller 
and  the  resulting  scaled  dose  is  within  ten  percent  of  the  dose  in  air  predicted 
by  CYLTRAN.  Figure  6.2  through  6.6  present  curves  of  the  dose  in  water 
scaled  to  air,  the  dose  in  liquid  nitrogen  scaled  to  air,  and  the  dose  in  air  as 
predicted  by  CYLTRAN. 

The  critical  energy  for  water,  air,  and  liquid  nitrogen  are  contained  in 
Table  1.  Since  the  experiment  and  the  calculations  were  conducted  at  100 
MeV,  close  to  the  critical  energies  of  the  three  media,  the  concept  of 
radiation  length  is  somewhat  nebulous.  Water  with  solid  angle  effects  can  be 
scaled  to  air  with  consistent  results  because  they  have  similar  Z.  It  is 
unlikely  the  results  in  water  would  scale  to  lead,  for  instance.  At  lower 
energies,  a  better  length  parameter  is  the  range  in  gm/cm2.  Figure  6.7  and 
Figure  6.8  contain  the  off-axis  dose  in  water.  Figure  6.9  and  Figure  6.10 
contain  the  results  in  liquid  nitrogen. 
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C.     COMPARISONS  WITH  OTHER  PAPERS 

In  his  CYLTRAN  predictions  for  liquid  nitrogen,  Cromar  included  a  0.2  cm 
layer  of  silicon  in  the  geometry  for  the  energy  deposition  [REF.  1 3 .  This 
special  geometry  required  individual  CYLTRAN  runs  for  each  distinct  length 
used  in  the  experiment,  requiring  a  great  deal  of  redundant  calculations  by 
the  computer.  To  avoid  this,  the  calulations  in  this  paper  were  all  done  in  a 
single  run  for  each  different  medium.  The  energy  deposition  was  recorded  in 
sub-zones  of  the  medium  instead  of  silicon.  To  ensure  the  validity  of  this 
approach,  the  off-axis  dose  in  liquid  nitrogen  predicted  by  this  method  is 
compared  to  the  experimental  results  reported  by  Cromar.  Figures  6.11 
through  6. 13  contain  the  two  sets  of  data  and  in  general  show  good  agreement 
between  the  predicted  and  observed  results. 

Cromar's  results  were  also  compared  to  those  reported  by  R.A  Lindgren 
[Ref.7]  which  looked  at  the  energy  depositon  in  air  due  to  50  MeV  electrons 
using  the  computer  code  ETRAN.  The  differences  between  the  normalized 
doses  were  significant  but  several  different  assumptions  were  made  in  the 
two  studies  and  the  difference  in  the  incident  energy  makes  a  direct 
comparison  difficult. 

An  additional  study  of  electron  transport  at  energies  up  to  1000  MeV  has 
been  reported  by  S.M.  Seltzer  and  M.J.  Berger  [Ref.8].  In  their  study,  the 
computer  code  ETRAN  was  used.  Since  the  TIGER  series  of  codes  all  draw  from 
the  same  source  for  the  Monte-Carlo  calculations  the  results  shoud  be  the 
same  as  those  presented  in  this  paper  generated  from  the  CYLTRAN  code. 
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Figure  6.14  shows  the  CYLTRAN  data  plotted  with  the  ETRAN  data  at  incident 
energies  of  60  MeV  and  125  MeV.  The  ETRAN  curves  were  derived  from 
replicate  data  of  Ref.  8  and  plotted  at  the  midpoints  of  the  histogram  bins. 
ETRAN  is  a  one-dimensional  code  and  for  comparison  it  was  necessary  to  sum 
the  ten  radial  sub-zones  used  by  CYLTRAN  to  obtain  a  single  value  of  energy 
deposition  at  a  given  penetration  depth.  The  CYLTRAN  zones  were  only 
0.2  g/cm2  in  thickness  so  an  additional  factor  of  five  was  included  to 
normalize  the  results  to  a  unit  depth.  Figure  6.14  shows  that  once  the 
geometries  are  equalized,  there  is  no  significant  differences  between  the  two 
computer  codes. 

D.    ENERGY  DEPENDENCE  OF  THE  RESPONSE 

An  additional  experimental  run  was  conducted  at  20  MeV  incident  energy 
however  the  results  were  not  returned  for  inclusion  in  this  paper.  For  future 
work,  it  is  recommended  that  the  Naval  Postgrtaduate  School  purchase  or 
obtain  the  use  of  its  own  TLD  reader  and  a  supply  of  reusable  TLDs.  This 
experiment  alone  used  80  TLDs  at  each  energy.  This  was  was  considered  the 
minimum  necessary  to  obtain  quick,  but  meaningful  results.  For  a  detailed 
study,  the  number  of  TLDs  could  easily  run  into  the  hundreds.  Having  to 
depend  on  a  third  party  to  process  this  number  of  TLDs  is  time  consuming  and 
imposses  a  huge  burden  of  the  processing  facility.  It  would  greatly  facilitate 
the  studies  if  the  TLDs  could  be  processed   at  the  test  facility. 


52 


E.     SUMMARY 

The  results  presented  indicate  that  the  CYLTRAN  computer  computations 
do  accurately  predict  the  experimental  results.  There  does  not  appear  to  be 
any  dlscreapncy  between  the  one-dimensional  computer  code  ETRAN  and  the 
two-dimensional  computer  code  CYLTRAN.  For  materials  of  similar  Z,  the 
extension  of  experimental  results  in  one  medium  to  another  appears  to  be 
valid. 
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TABLE  6.1 
PREDICTIONS  FOR  AIR  BASED  ON  RADS/COULOMB  IN  WATER 

Factor=(30700/37T2 


Rad/Lth 


Air  (Pred)  One-Slgma    Air  (Cyltran)      One-Slgma     Ratio 

Rads/Coulomb     Error  Rads/Coulomb     Error 


0.00 


7.9245E+04 

2 

9.2801E+04 

3 

1.17 

8.1257E-01 

56 

3.1579E+00 

91 

3.89 

7.5299E-01 

46 

3.0577E-01 

42 

0.41 

4.6982E-01 

48 

2.2514E+00 

62 

4.79 

9.9008E-01 

30 

I.0680E+00 

42 

1.08 

8.3157E-01 

28 

3.I993E-01 

34 

0.38 

1.2156E+00 

27 

3.8732E-01 

43 

0.32 

5.5566E-01 

20 

4.8336E-01 

36 

0.87 

4.0331E-01 

24 

2.4460E-01 

47 

0.61 

2.1545E-01 

27 

1.8029E-01 

51 

0.84 

0.25 


8.2900E+04 

3 

9.5737E+04 

3 

1.15 

4.0395E+03 

3 

3.1162E+03 

6 

0.77 

3.2665E+02 

7 

3.0579E+02 

9 

0.94 

8.0173E+01 

13 

8.9421E+01 

18 

1.12 

3.8481E+01 

16 

3.2072E+01 

28 

0.83 

1.2797E+01 

25 

4.1353E+00 

99 

0.32 

5.5583E+00 

26 

8.731 5E+00 

40 

1.57 

3.4629E+00 

27 

4.5326E+00 

•  40 

1.31 

3.1313E+00 

27 

3.2345E+00 

55 

1.03 

2.8272E+00 

52 

3.2645E+00 

61 

1.15 
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TABLE  6.1    (cont'd) 


0.50 


2.6881E+04 

3 

2.9248E+04 

3 

1.09 

1.2841E+04 

2 

1.3804E+04 

2 

1.07 

4.1262E+03 

3 

4.2394E+03 

4 

1.03 

1.1693E+03 

5 

1.1804E+03 

6 

1.01 

3.7070E+02 

4 

3.4757E+02 

7 

0.94 

1.0736E+02 

13 

1.1207E+02 

17 

1.04 

5.5645E+01 

8 

5.5706E+01 

11 

1.00 

3.2064E+01 

11 

3.5135E+01 

15 

1.10 

9.5732E+00 

28 

1.2687E+01 

17 

1.33 

8.5832E+00 

15 

8.8687E+00 

29 

1.03 

00 
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2.3485E+03 

6 

2.2181E+03 

7 

0.94 

1.7406E+03 

6 

1.6862E+03 

6 

0.97 

1.3403E+03 

3 

1.1146E+03 

6 

0.83 

9.5150E+02 

3 

7.4772E+02 

5 

0.79 

6.8564E+02 

4 

5.8765E+02 

6 

0.86 

4.8474E+02 

7 

3.5625E+02 

4 

0.73 

3.4158E+02 

2 

2.4641E+02 

9 

0.72 

2.2934E+02 

7 

1.7884E+02 

5 

0.78 

1 . 5626E+02 

3 

1.1710E+02 

8 

0.75 

8.9038E+01 

12 

8.3028E+01 

10 

0.93 

5.3640E+02 

10 

5.5813E+02 

16 

1.04 

3.7059E+02 

7 

4.4524E+02 

5 

1.20 

3.0896E+02 

8 

3.0862E+02 

5 

1.00 

1.7579E+02 

12 

1.8005E+02 

7 

1.02 

1.1369E+02 

10 

1.2701E+02 

14 

1.12 

1.0442E+02 

7 

7.9514E+01 

10 

0.76 

7.0638E+01 

13 

6.0580E+01 

9 

0.86 

5.6338E+01 

12 

4.2238E+01 

15 

0.75 

4.2163E+01 

12 

3.0173E+01 

8 

0.72 

2.4319E+01 

13 

2.5074E+01 

18 

1.03 
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APPENDIX  A 
DESCRIPTION  OF  CYLTRAN  AND  OUTPUT 

The  computer  code  CYLTRAN  is  a  member  of  the  Integrated  TIGER  Series  of 
Coupled  Electron/Photon  Monte  Carlo  Transport  Codes  developed  at  Sandia 
National  Laboratories  by  J. A.  Halbleib  and  W.H.  Vandevender.  CYLTRAN 
employs  a  fully  three-dimensional  description  of  particle  trajectories  within 
an  axisymmetric  cylindrical  material  geometry  and  finds  application  in 
problems  involving  electron  and  photon  beam  sources. [Ref. 6:  p. 6] 

CYLTRAN  is  based  primarily  on  the  ETRAN  model  which  combines 
microscopic  photon  transport  with  a  macroscopic  random  walk  for  electron 
transport.  The  code  describes  the  generation  and  transport  of  the 
electron/photon  shower  from  the  incident  energy  down  to  a  specified  cut-off 
energy  for  photons  and  electrons.  The  code  can  handle  up  to  ten  different 
materials,  with  the  only  requirement  being  that  the  material  geometry  be 
cylindrical ly  symmetric. [Ref. 6:  p. 70] 

The  user  may  specify  incident  particles  to  be  either  electrons  or  photons 
which  may  be  considered  mono-energetic  or  have  a  specified  energy 
spectrum.  The  source  reference  direction  can  be  monodirectional , 
isotropic,  or  have  a  specified  divergence  angle. 
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The  default  output  of  CYLTRAN  consists  of: 

a)  Energy  and  number  escape  fractions  (leakage)  for  electrons, 
unscattered  photons  and  scattered  photons. 

b)  Charge  and  energy  deposition  profiles. 

c)  An  explicit  statement  of  energy  conservation. 

In  addition  to  the  default  output,  a  number  of  optional  outputs  may  be 
selected  through  the  use  of  the  appropriate  keywords.     These  are: 

a)  Escape  fractions  that  are  differantial   in  energy  for  both  electrons 
and  scattered  photons.     . 

b)  Escape  fractions  that  are  differential  in  angle  for  both  electrons  and 
scattered  photons. 

c)  Coupled  energy  and  angular  distributions  of  escaping  electrons  and 
scattered  photons. 

d)  Volume-averaged  energy  distributions  of  electron  and  photon  scalar 
fluxes  for  selected  regions  of  the  problem  geometry. 

e)  Psuedo-pulse-height      distributions      for      selected   regions   of  the 
problem  geonetry —  for  example,   those  regions     corresponding     to 
active  detector  elements. [Ref. 6:  pp.  19,34] 

Two  basic  steps  are  required  for  solving  any  given  problem: 

a)  Generate  the  cross  sections  by  running  the  cross  section  code. 

b)  Run  the  Monte  Carlo  code. 

Table  Al  is  a  sample  input  used  for  the  generations  of  cross-sections  of 
water,  standard  air  (80%  nitogen,  20%  oxygen),  liquid  nitrogen,  and  silicon. 
Table  A2  is  a  sample  input  for  the  proper  execution  of  CYLTRAN  once  the 
cross-sections  have  been  generated.  CYLTRAN  output  for  runs  of  20,000 
incident  electron  in  water,  air,  and  liquid  nitr.ogen  is  include  in  Tables  A3 
thru  Table  A25. 
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TABLE  A1 
SAMPLE  INPUT  FOR  THE  GENERATION  OF  CROSS  SECTIONS 

1  MATERIAL  H  .  I  11  0  .  889  DENSITY  1.0 

2  MATERIAL  N  .80  0  .20  GAS  DENSITY  .001226 

3  MATERIAL  N  DENSITY  0.81 

4  MATERIAL  SI 

5  TITLE    100  MeV  CROSS  SECTIONS  FOR  H20  AIR  N  SI 

6  ENERGY  100.0 


'2 


TABLE  A2 


SAMPLE  INPUT  FOR  THE  EXECUTION  OF  CYLTRAN 


1  ECHO  1 

2  TITLE 

3  ...  100  MeV  ENERGY  DEPOSITION  IN  WATER 

4  ****** ****************S0URCE ******************** ******* 

6  ENERGY  100.0 

7  *  DEFAULT  SOURCE  PHASE  SPACE  PARAMETERS 

8  POSITION  0.0  0.0  0  0.0 

9  RADIUS  0.4 

10  *  MONODIRECTIONAL  Z-AXIS 

11  DIRECTION  0.0  0.0  0.0 

1 2  ************ ********** OPT  IONS************************** 

13  HISTORIES  20000 

14  BATCHES  10 

15  CUTOFFS  0.05  0.001 

16  DUMP 

17  *********************0UTPUT  OPTIONS******************** 
ELECTRON-ESCAPE 

NBINE  2 

NBINT  4 
PHOTON-ESCAPE 

NBINE  2 

NBINT  4 
ELECTRON-FLUX 
NBINE  2 

26  PHOTON  -FLUX  1 

27  NBINE  2 

28  **************** ****GEOMETRY* ************************* 


18 
19 
20 
21 
22 
23 
24 
25 


10 


0 


29 

GEOMETRY  11 

30 

*ZL 

ZR 

Rl 

RO    MAT 

31 

0.0 

0.2 

0.0 

10.0  1 

1   1   10 

32 

0.2 

9.2 

0.0 

10.0  1 

33 

9.2 

9.4 

0.0 

10.0  1 

1   1   10 

34 

9.4 

18.5 

0.0 

10.0  1 

35 

18.5 

18.7 

0.0 

10.0  1 

1   1   10 

36 

18.7 

37.0 

0.0 

10.0  1 

37 

37.0 

37.2 

0.0 

10.0  1 

1   1   10 

38 

37.2 

55.5 

0.0 

10.0  1 

39 

55.5 

55.7 

0.0 

10.0  1 

1   1   10 

40 

55.7 

74.0 

0.0 

10.0  1 

NZ   NR   ECUT   PTCZ 
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APPENDIX  F 


EQUATIONS  USED  FOR  CURVE  FITTING 


TABLE  F1  :  100  MEV  ELECTRONS  IN  WATER 


0.0  cm 


0.5i  x  i  1.5 
1.5<  x   i  10 


y= 1 . 7 1  x  1 0 ' ' 

y=7. 516       x  105 


-11.49x 
[-0.106IX 


9.2  cm 


0.5i  x  i  1.5 
1.5<  Xi  2.5 

2.5<  x  i  3.5 
3.5<  x  i  6.5 
6.5<  x  i 


7.5<  x  i 


7.5 

10 


y=2.59 

X 

y=1.211 

X 

y=7.514 

X 

y= 1.359 

X 

y=8.291 

X 

y=5. 1014 

X 

10 


11 
11 


10* 
10< 
10' 

10< 


,-3.022x 
[-2.515X 

[-1.405X 
;-0.9046x 
;~0.4732x 
[-O.IOHx 


18.5  cm 


0.5i  x  i  1.5 
1.5<  x  i  2.5 
2 .  5<  x  i  3 .  5 
3. 5<  x  i  4.5 
4.5<  x  i  5.5 
5.5<  x  i  10 


y=2.6675 
y=4.8532 
y=6.6444 
y=4.4876 
y=6.7422 
y=2.2332 


xlO10 

e-0.7388x 

xlO10 

e-1. 1352x 

xlO10 

e-l. 1352x 

xlO10 

e-1 . 1488x 

xio'0 

e-1.2393x 

xlO9 

e-0.629x 

37.0  cm 


0.5i  xi  10 


y=1  .848-4.871  x  108x  +4.807  x  107x2 
1.748  x  106x3 


55.5  cm 


0.5i  x  ilO 


y=4.26         xlO8  e-°-3299x 


74.0 


0.5i  x  i  10 


y=2.074        x  10 


8 


1-0.282x 
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TABLE  F2  :  100  MEV  ELECTRONS  IN  AIR 


0.0  m 


4iXi  12 

12<  x   i  80 


y= 1.587 
y=5.972 


x 
x 


>-1.24x 

[-0.0446X 


76.8  m 


53.5  m 


4i  x  i  12 

y=5.297 

x  105 

e-0.4126x 

12<  Xi  20 

y=1.017 

x  105 

e-0.2792x 

20<  x  i  28 

y=6.614 

x  103 

e-0.1481x 

28<  x  i  36 

y=3.236 

x  103 

e-0.1253x 

36<  x  i  54 

y=5.50 

x  102 

e-0.07838x 

54<  x  i  80 

y=215.3-8 

.25x+0. 

1066x2-.000458x3 

4iXi  12 

y=4.262 

x  104 

e-0.0905x 

12<  xi  20 

y=8.1107 

xlO4 

e-0.1422x 

20<  xi  28 

y= 1.0363 

x  105 

e-0. 1540x 

28<  x  i  36 

y=8.5209 

x  104 

e-0. 1473x 

36<  x  i  44 

y=5.6628 

x  104 

e-0.1364x 

44<  x  i  80 

y=4. 1004 

x  103 

e-0.0789x 

307.0  m 


4i  x  i  80 


y=2.802       x  10- 


e-0.04468x 


460.5  m 


4i  x  i80 


y=6.90         x  10' 


>-0.04427x 
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0.0  cm 


TABLE  F3  :  100  MeV  Electrons  in  Liquid  Nitrogen 


0.5i  x  i  1.5 
1.5<  x  i  10 


y=6. 166       x  1012  e"n-49x 

y=4.198       xlO6  e"°-1061x 


1 1.9  cm 


0.5ixi2.5 
2.5<  xi  3.5 
3.5<  x  i  4.5 
4.5<  xi  10 


y=  1.593       x  1011  q-3.022* 

y=2.162       xlO10 
y=5. 198       x  109 


y=5.525       x  10 


■2.515x 
^-1.405x 


8  a-0.9046x 


23.8  cm 


0.5i  x  i  1 .5 
1.5<  xi  2.5 
2.5<  x  i  6.5 
6.5<  xi  10 


10  a-1.1352x 


y= 1.369  xlO10  e~°-7388x 

y=2.397  xio" 

y=3.873  x  10 

y- 1.853  xlO10  e-'-1488x 


10 


1352x 


47.7  cm 


0.5i  xi  10 


y-1.17         x  109  e-°-258x 


71.5  cm 


0.5i  x  ilO 


y=2.61  xlO8  e-°-28,7x 


95.4 


0.5i  x  i  10 


y-1.221       xlO8  e-°-2319x 
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TABLE  F4  :  20  MeV  Electrons  in  Water 


0.0  cm 


0.5i  x  i 
1 .5<  x  i 


.5 
0 


y=5.4289 
y= 1 . 9 1 99 


x  10p 
x  10' 


12 


-9.1981X 
!-0.4918x 


9.2  cm 


0.5iXi  1.5 
1.5<  xi  2.5 
2.5<  x  i  3.5 
3.5<  x  i  4.5 
4.5<  x  i  6.5 
6.5<  x  i  7.5 
7.5<  x  i  8.5 
8.5<  x  i  10 


y=4.1483 
y=S.3790 
y=2.1945 
y=1 .7157 

y= 1.6720 
y=6.9153 
y=5.8735 
y=5.0129 


xlO9 
xlO9 


10 
11 


x 

xlO10 
x  108 


10 
I0< 


3 


-0.3287X 
-0.7938X 
1825x 
-1.77x 
-1.262x 
-0.7654x 
-0.7436x 
-0.1832x 


18.5  cm 


0.5iXi  10 


y=6.5584  x  107-  2.  1939  x  107x  + 


2.6384  x  106x2  - 


0869  x  105x3 


37.0  cm 


0.5iXi  10 


y=7.572       x  10( 


■0.209x 


55.5  cm 


0.5i  x  ilO 


y=8.286       x  10( 


,-0.3157x 
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APPENDIX  6 
EXPERIMENTAL  RESULTS 

Table  Gl  contains  the  results  of  the  100  MeV  experimental  run.  The 
uncertainties  involved  include  the  TLD  position  on  the  stretcher  (±  0.5  cm), 
the  strecher  location  with  respect  to  the  central  axis  of  the  test  tank  (±  1 .0 
cm),  and  the  skewness  of  the  beam  with  respect  the  center  line  of  the  tank 
(±2.0  degrees).  Uncertainty  in  the  charge  delivered  is  five  percent  and 
uncertainty  in  the  normalized  dose  is  ten  percent. 
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TABLE  G1 
Experimental  Run  #1:  100  Mev  Electrons  in  Water 

Test  day -28  JUNE  1985 

Test  Medium  -  Water 

Beam  Diameter  -  0.5  cm 

(0.4  cm  above  the  horizontal  axis) 

SEM  Efficiency-  0.028 

Capacitor-  IxlO"6  Farad 

TLD  Type  -  CaF2 


TLD 

z-  axis 

x-axis 

Rads 

Coulomb 

Rads/Coul 

Number 

(cm) 

(cm) 

1 

0.0 

2.0 

35.6702 

1.03E-06 

3.46E+07 

2 

0.0 

4.0 

15.8316 

1.03E-06 

1.54E+07 

3 

0.0 

6.0 

64.0929 

1.36E-05 

4.71E+06 

4 

0.0 

8.0 

49.6843 

1.36E-05 

3.65E+06 

5 

9.2 

1.0 

2348.76 

1.03E-06 

2.28E+09 

6 

9.2 

2.0 

631.758 

1.03E-06 

6.13E+08 

7 

9.2 

3.0 

183.113 

1.03E-06 

1.78E+08 

8 

9.2 

4.0 

78.6993 

1.03E-06 

7.64E+07 

9 

9.2 

5.0 

42.9735 

1.03E-06 

4. 17E+07 

10 

9.2 

6.0 

22.6433 

1.03E-06 

2.20E+07 

11 

9.2 

7.0 

18.3424 

1.03E-06 

1.78E+07 

12 

9.2 

8.0 

111.069 

1.36E-05 

8. 17E+06 

13 

9.2 

9.0 

111.267 

1.36E-05 

8. 18E+06 
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TABLE  G1  (cont'd) 


TLD 

z-  axis 

x-axis 

Rads 

Coulomb 

Rads/Coul 

Number 

(cm) 

(cm) 

14 

18.5 

1.0 

12188.1 

1.03E-06 

1.18E+10 

15 

18.5 

2.0 

4085.73 

1.03E-06 

3.97E+09 

16 

18.5 

3.0 

2407.99 

1.03E-06 

2.34E+09 

17 

18.5 

4.0 

516.688 

1.03E-06 

5.02E+08 

18 

18.5 

5.0 

297.201 

1.03E-06 

2.89E+08 

19 

18.5 

6.0 

97.8459 

1.03E-06 

9.50E+07 

20 

18.5 

7.0 

453.528 

1.36E-05 

3.33E+07 

21 

18.5 

8.0 

340.428 

1.36E-05 

2.50E+07 

22 

18.5 

9.0 

194.366 

1.36E-05 

1.43E+07 

23 

18.5 

10.0 

156.074 

1.36E-05 

1.15E+07 

24 

37.0 

1.0 

1282.73 

1.03E-06 

1.25E+09 

25 

37.0 

2.0 

1783.29 

1.03E-06 

1.73E+09 

26 

37.0 

4.0 

3612.04 

1.36E-05 

2.66E+08 

27 

37.0 

6.0 

5309.53 

1.36E-05 

3.90E+08 

28 

37.0 

8.0 

4737.13 

1.02E-04 

4.64E+07 

29 

37.0 

10.0 

7204.42 

1.02E-04 

7.06E+07 

30 

55.5 

1.0 

5743.79 

1.36E-05 

4.22E+08 

31 

55.5 

2.0 

1189.37 

1.36E-05 

8.75E+07 

32 

55.5 

4.0 

13359.7 

1.02E-04 

1.31E+08 

33 

55.5 

6.0 

3039.64 

1.02E-04 

2.98E+07 

34 

55.5 

8.0 

4717.41 

1.02E-04 

4.62E+07 

35 

55.5 

10.0 

1047.85 

1.02E-04 

1.03E+07 

36 

74.0 

2.0 

2625.15 

1.36E-05 

1.93E+08 

37 

74.0 

4.0 

2625.15 

1.02E-04 

2.57E+07 

38 

74.0 

6.0 

7757.14 

1.02E-04 

7.61E+07 

39 

74.0 

8.0 

1181.88 

1.02E-04 

1. 16E+07 
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